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The flow properties of a continuum model for an active nematic is studied and compared with 
recent experiments on suspensions of microtubule bundles and molecular motors. The velocity cor- 
relation length is found to be independent of the strength of the activity while the characteristic 
velocity scale increases monotonically as the activity is increased, both in agreement with the ex- 
perimental observations. We argue that the length scale is set by the density of topological defects. 



Active systems that produce their own energy, such 
as bacterial suspensions and living cells, are proving rich 
research areas [1, 2]. Apart from their obvious relevance 
to a quantitative understanding of biological functions, 
active systems naturally operate out of thermodynamic 
equilibrium and hence provide a testing ground for theo- 
ries of non-equilibrium statistical physics [3] . Moreover it 
may be feasible to harness the chemical energy of molecu- 
lar motors or the motion of microswimmcrs to drive novel 
micromachincs [4-6]. 

Hydrodynamic instabilities are inherent to active flu- 
ids [7-10], and active suspensions spontaneously gener- 
ate complex flow patterns at length scales much larger 
than the actual constituents [11-17]. These are char- 
acterised by strong variations in vorticity (see Fig. 2). 
Such turbulent-like patterns have been observed in ex- 
periments on mixtures of actin or microtubules and mo- 
tor proteins designed to isolate the important ingredients 
leading to cellular motility [6, 18-20]. Very similar struc- 
tured flows have been observed in two-dimensional layers 
of epithelial cells [21] and, at larger length scales, in dense 
suspensions of swimming bacteria [11, 12, 16]. The fea- 
tures of the flow have been reproduced qualitatively in 
simulations of driven rods [22, 23] and microswimmers 
[13, 17, 24] and in continuum theories of active nematics 
[15, 16, 25]. 

The similarity in the hydrodynamic behaviour of the 
different active experimental and model systems is very 
appealing, but there are still large gaps in our under- 
standing of the physical mechanisms driving active tur- 
bulence. Moreover the extent of and reasons behind any 
'universal' behaviour remains unclear [16]. To help ad- 
dress these questions in this paper we measure the root 
mean square (RMS) velocity and velocity- velocity corre- 
lation function in a continuum model of an active sus- 
pension, as a function of the activity, and compare it 
to recent experiments on extensile microtubule bundles 
driven by the motor protein kincsin [6] . This comparison 
is summarised in Fig. 1. 

Figure 1(a) shows that the normalised velocity- velocity 
correlation functions collapse onto a single curve, inde- 
pendent of the strength of the activity over a wide range 
of activities. The same scaling behaviour is seen in the 
experiments (see Fig 2(c) of Rcf. [6]). This indicates 
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FIG. 1. The velocity field generated in a two dimensional ac- 
tive nematic suspension is characterised by (a) the velocity- 
velocity correlation function and (b) the RMS velocity. In 
Fig. 1(a), normalised velocity- velocity correlations are plotted 
as a function of distance, R for different values of the activ- 
ity. Unnormalised functions are shown in the inset (i). The 
collapse of the normalised correlation functions onto a single 
curve indicates insignificant correlation between the charac- 
teristic flow length scales and activity. This is also illustrated 
in inset (ii) where a characteristic length, defined as the value 
of R where the velocity- velocity correlation function has fallen 
to 0.8 of its maximum value, is plotted as a function of ac- 
tivity, C Fig. 1(b) shows a comparison between experimental 
[6] and simulation results for the variation of the RMS veloc- 
ity with activity. Mapping of simulation to physical units for 
velocity is explained in the text. 



that there is a length scale in the problem which is set 
by the structure of the underlying fluid, independent of 
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the strength of the activity. Figure 1(b) compares the 
RMS velocity in the simulations and experiments where, 
to make the comparison, we assume a logarithmic depen- 
dence of the activity on ATP concentration. Two distinct 
regions are observed, a gentle increase in the velocity with 
activity at low activities, followed by a sharper increase 
at higher activities. (The simulations then become un- 
stable and we do not see the saturation observed in the 
experiments.) 

We will relate this behaviour to features of the velocity 
field, which we interpret in terms of velocity jets driven 
by defects and distortions of the director field [10]. The 
length scale of velocity correlations is set by the density 
of defects which depends strongly on the value of the 
elastic constant but only very weakly on the strength of 
the activity. The magnitude of the velocity, however, de- 
pends on the flow induced by the defects, which increases 
linearly with activity. At very low activities, individual 
defects are replaced by lines of kinks in the director field, 
which lead to smaller velocity magnitudes. We show that 
similar features are seen for both extensile and contractile 
active nematics, despite the coupling between the defects 
and the flow field being different in the two cases. 

First we describe the equations of motion, those corre- 
sponding to an active nematic, that we use to model the 
active suspension. We represent the director field using 
the alignment tensor Q. Evolution of Q along with the 
momentum pu is given by [26, 27], 

(9t + U fc 0fc)Qy-5y=rJ?y, (1) 

p(d t + u k d k )ui = djUij. (2) 

The strain rate tensor, Eij = (diUj + djUi)/2 and the 
vorticity tensor, Ojj = (djUi — diUj)/2 describe the gen- 
eralised advection term 

Sij =(XE lk + n ik )(Q kj + S kj /3) + (Q lk + <W3) 
(AE kj - Sl kj ) - 2A(Q 4 , + 6 l3 /3)(Q kl d k ui), 

where A is the alignment parameter. We choose A = 0.7 
corresponding to tumbling rods. Rotational diffusivity is 
denoted by T and the molecular field = —5J-/8Qij + 
(Sij ■ /3)Tt(SJ- / 5Q k i) is determined from the free energy: 

T = K(d k Q tJ ) 2 /2 

+ AQijQji/2 + BQ i:j Q jk Q ki /'S + C(Q tJ Q 3l ) 2 /A. (3) 

Here K is the elastic constant, A, B and C are material 
constants. The total stress generating the hydrodynam- 
ics has 3 parts; (i) the viscous stress, U^ scous = 2p,Eij, 
(ii) the passive stress, 

mr me = ~ P <% + 2 KQij + 5ij/S){QkiH lk ) 

- \H lk (Q k3 + S kj /3) - \(Q ik + 5 ik /3)H kj 

— diQ k i h Qi k H k j — H lk Q kj (4) 

SdjQi k 



and the (iii) active stress, ]J99 tlve = —QQij introduced 
in [7]. Here ( is called the activity coefficient. Note 
that P = Pop — K(d k Qij) 2 /2 is a good approximation 
to the modified hydrostatic stress [28]. We solve this 
coupled system of governing equations (Eqs. (1) and (2)) 
using a hybrid lattice Boltzmann method on a D3Q15 
lattice and we refer the reader to [29-31] for the details 
of the equations and solution method. Simulations are 
performed on a 400 x 400 domain and the parameters used 
are T = 0.34, A = 0.0, B = -0.3, C = 0.3, K = 0.02, 
Pq = 0.25 in lattice units, unless mentioned otherwise. 

It is well known that active liquid crystal equations 
predict hydrodynamic instabilities [7-9] leading to irreg- 
ular, turbulent-like flow patterns. The left column of 
Fig. 2 depicts a typical flow pattern observed for an ex- 
tensile suspension. The top figure is on the scale of ne- 
matic correlation lengths, the middle frame corresponds 
to intermediate lengths at which individual vortices are 
visible, and bottom frame to hydrodynamic length scales. 
The region within the white border in each figure corre- 
sponds to the figure above. In each figure, colour shading, 
continuous lines and dashes correspond to the vorticity, 
stream lines and director field respectively. 

Consider first Fig. 2(a). This shows a topological de- 
fect of charge +1/2 in the director field. The director field 
consequently has a strong bend deformation which leads 
to an active stress that generates a flow [10] along the 
axis of symmetry of the defect, indicated by the stream- 
lines in the figure. The vorticity field helps to locate the 
source of the jet, as a transition from blue, represent- 
ing clockwise vorticity, to red, representing anticlockwise 
vorticity. 

Moving to larger length scales, a large number of de- 
fects, and hence jets, may be identified, most easily by 
sharp changes in colour of the vorticity field. The posi- 
tion and orientation of the defects are random, leading 
to jets that interact to produce complex velocity fields. 
For example, Fig. 2(b) shows four +1/2 defects that lead 
to a circulating velocity field. At the larger length scale 
of Fig. 2(c) several such circulating flows may be seen 
each of which can be identified with two or more bend 
defects that produce jets moving in different directions. 
Defects oriented in the same direction generate paral- 
lel flows and hence longer and stronger jets. Defects of 
charge — 1/2 can also be identified in the field. These may 
be identified as a combination of three bend or splay de- 
formations oriented such that the jets produced by each 
tend to cancel out. Therefore such defects act as effective 
stagnation points, mostly altering the flow directions of 
jets and hence contributing to the circulating flows. 

The right hand column of Fig. 2 presents similar data 
for a contractile suspension. Again the flow is driven 
by variations in the director field. However, in the ex- 
tensile case these are predominantly bend, distortions, 
whereas in the contractile case they are predominantly 
splay. This corresponds to the primary instabilities in 
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FIG. 2. Flow and director fields at three different length scales 
in active nematic suspensions. The length scale increases from 
top to bottom, and each frame occupies the white square in 
the frame below. Figures in the left column (right column) 
correspond to extensile (contractile) stress in the suspension. 
Colour shading, continuous lines, and dashes represent vor- 
ticity, stream lines and the director field respectively. Defects 
and deformations in the director field generate fluid jets as 
illustrated in Fig. (a), (d). (Note that a defect source can 
be distinguished by a sharp colour change in the vorticity. 
Blue and red colours correspond to clockwise and anticlock- 
wise vorticity respectively.) Two or more jets lead to circulat- 
ing flows (Figs, (b), (e)) resulting in the complex flow patterns 
of Figs, (c), (f). 

the active nematic director field which are bend for a 
contractile, tumbling material and splay for an extensile 
one [9, 10]. The velocity fields are ~ 50 — 70% lower 
in the contractile suspensions for the same value of the 
activity. Indeed, the critical activity needed to trigger 
spontaneous active flow in one dimension is known to be 
lower for an extensile than for a contractile system [10]. 
The circulating flows can be easily distinguished from 



the vortices of high Reynolds number (Re) turbulence. 
At high Re an energy input at the largest scale causes a 
laminar-turbulent transition generating vortices through 
the energy cascading mechanism [32]. In contrast, in ac- 
tive fluids the energy input at the smallest scale is trans- 
ferred to the largest scale through the formation of jets 
and hence vortices. 

Having characterised the velocity field in terms of de- 
fects and jets we can motivate why the decay of the nor- 
malised velocity-velocity correlation function is indepen- 
dent of the activity. This follows because the defect den- 
sity controls the spatial structure of the flow field, and 
hence the dependence of the velocity correlations on dis- 
tance. Indeed, we find that the number of defects [33] 
in the simulation box depends only very weakly on the 
strength of the activity for fixed values of the elastic con- 
stant, as shown in Fig. 3. Although the number of +1/2 
defects is shown, the number of —1/2 defects is the same, 
further indirect evidence that the liquid crystal free en- 
ergy, rather than the different interaction of the ±1/2 
defects with the flow field, dominates in controlling their 
number density. Moreover, both extensile and contractile 
suspensions have same number of defects. 
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FIG. 3. Total number of point 1/2 defects plotted as a func- 
tion of (i) activity for extensile and contractile suspensions, 
(ii) elasticity for extensile suspensions. The number of point 
defects decreases strongly with increase in the elasticity of the 
material, while there is only a weak dependence on activity. 

To further check that the defect density controls the 
length scale of the flow field we also measured the num- 
ber of point defects within the domain as a function 
of the elastic constant. Figure 3 shows that the num- 
ber of defects increases sharply as the elastic constant 
is decreased. The normalised velocity-velocity correla- 
tion functions, plotted in Fig. 4, do not scale for different 
elasticity, but decay more rapidly with decreasing clastic 
constant corresponding to the smaller spacing between 
the topological defects. 

Although the spatial scale of the flow is indepen- 
dent of activity, higher activity drives stronger jets 
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FIG. 4. Normalised velocity-velocity correlation functions as 
a function of distance for different values of the elastic con- 
stant. Non-normalised functions are plotted in the inset. A 
monotonic decrease in the length scale is seen as the elasticity 
of the material is reduced. 



[10] and hence the RMS velocity of the flow field in- 
creases, approximately linearly with increasing activity. 
At very small activities Fig. 1 indicates that there is 
a regime where the flow velocity increases significantly 
more slowly. Similar behaviour is seen in the micotubule- 
kinesin experiments, although no comment is made in the 
paper [6]. Moreover, if concentration of swimmers may 
be considered as a measure of activity, a similar crossover 
was observed in particle simulations of extensile swim- 
mers [17] where it was interpreted as a transition from 
a weakly correlated to strongly correlated dynamics of 
swimmers. In our simulations the low activity regime 
corresponds to the existence of lines of kinks in the di- 
rector field, very similar to those observed in stable one- 
dimensional flows [29] . These structures are illustrated in 
Fig. 5 for an extensile and contractile suspension. The re- 
sults suggest that such low values of the activity are suffi- 
cient to create, but not to further break up, the quasi- ID 
kinked director and flow configurations. 

Our observed results could be put in perspective us- 
ing a simplistic scaling argument. The dynamics of the 
director field and the velocity field are interconnected 
through a feedback loop: nematic order creates an ac- 
tive stress that leads to flow, and the flow field in turn 
advects the nematic order. This leads to the instability 
that manifests itself through the proliferation of defects, 
leading to a dynamic structure that resembles the high 
temperature phase of a system with a defect-mediated 
phase transition, such as the Kosterlitz-Thouless transi- 
tion [34]. In such a regime, the freely fluctuating defects 
are expected to lead to Debye screening for the nematic 
order parameter, with a screening length I that depends 
on the average density of the defects n via I ~ n~ x / 2 . 
Balancing the rate of energy input via the active stress 



FIG. 5. Patterns formed at low activity (£ = 0.001) in (a) 
extensile and (b) contractile suspensions. Colour shading in- 
dicates the vorticity field. The director field is superimposed 
on one quarter of the plot, and streamlines are shown else- 
where. At these smaller activities (regime I in Fig. 1(b)) lines 
between regions of different vorticity mirror lines of kinks in 
the director field. 

and the rate of viscous dissipation over the domain sur- 
rounding every single defect, we find v ~ ££Q//j,, which 
confirms our finding for the dependence of the velocity 
scale on the activity. Moreover, the scaling suggests that 
(vv) depends on (QQ), and in particular, the correlation 
length for the velocity decay is controlled by the nematic 
screening length. 

A quantitative comparison of experiment and simula- 
tion needs us to choose a length scale L and a time scale 
T . L was measured from the decay of the (vv) correla- 
tions. T was estimated by choosing the activity at which 
the transition in slope occurs in Fig. 1(b) as a reference 
value and, for this value, using passive scalar tracking to 
obtain the crossover time-scale from the ballistic to the 
diffusive regime. This gives T = 200s and L = 200fim for 
the experiments (see Fig. 2 in [6]) and 100 and 10 4 lattice 
Boltzmann units for the simulations, and hence charac- 
teristic velocities of lfim/s and 0.01 LB units. These val- 
ues are used to normalise the RMS velocities in Fig. 1(b) 
giving very similar numerical values. 

To summarise, we have shown striking similarities be- 
tween the properties of a suspension of microtubule bun- 
dles and kinesin and the continuum theory of an active 
nematic, allowing us to make an estimate of the length 
and time scales that map the theory to the physical sys- 
tem. Thus, this study abridges the existing gap between 
theory and experiments of active nematics. In particular 
the decay of the velocity- velocity correlations with dis- 
tance is independent of activity, and the RMS velocity 
shows two regimes, increasing at first slowly, but then 
faster, with the activity. In the active nematic the flows 
are generated by defects and distortions in the director 
field, and the number of defects which is essentially inde- 
pendent of the activity, sets the length scale of the flow 
field. The low activity regime corresponds to one di- 
mensional lines of distortion in the director field, which 
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break up in to defects at higher activities. Both 1/2 
and —1/2 defects have been observed experimentally [6] 
and it would be interesting to relate their density to the 
length scale set by the velocity correlations. We also hope 
to motivate studies of active turbulence in other systems 
to look for similar properties of the flow field. 
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mental data in Fig. 1, and to M. Blow, J. Dunkel and M. 
Ravnik for helpful advice. We acknowledge funding from 
the ERC Advanced Grant MiCE. 
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